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I. INTRODUCTION 

In this article we report results from studies of the 
7*7 — > P transition form factors, where P is a pseu- 
doscalar meson. In our previous works [l|, Q the two- 
photon-fusion reaction 

e+e" ^ e+e"P, (1) 

illustrated by Fig. [TJ was used to measure the 7r° and 
?7c transition form factors. Here this technique is ap- 
plied to study the ?/ and ?/ form factors. The transition 
form factor describes the effect of the strong interaction 
on the 7*7* — )- P transition. It is a function, F{q{,q^), 



*Now at Temple University, Philadelphia, Pennsylvania 19122, 
USA 
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of the photon virtualities qf. We measure the differen- 
tial cross sections for the processes e+e^ — > e^e'-rf-'^ 
in the single tag mode where one of the outgoing elec- 
trons^ (tagged) is detected while the other (untagged) is 
scattered at a small angle. The tagged electron emits a 
highly off-shell photon with the momentum transfer = 
— = ip — p')^ , where p and p' arc the four-momenta 
of the initial and final electrons. The momentum trans- 
fer to the untagged electron (g|) is near zero. The form 
factor extracted from the single tag experiment is a func- 
tion of one of the q^'s: F{Q^) = P(-Q^O). To relate 
the differential cross section da{e'^e~ — >■ e~^e~P)/dQ'^ to 
the transition form factor we use formulae equivalent to 
those for the e+e~ —J- e"''e~7r° cross section in Eqs. (2.1) 
and (4.5) of Ref. 

At large momentum transfer, perturbative QCD pre- 
dicts that the transition form factor can be represented 
as a convolution of a calculable hard-scattering ampli- 
tude for 77* — ^ qq with a nonperturbative meson dis- 



^ Unless otherwise specified, we use the term "electron" for either 
an electron or a positron. 
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FIG. 1: The diagram for the e+e — >■ e e P two-photon 
production process, where P is a pseudoscalar meson. 

tributioii amplitude (DA) (j)p{x,Q'^) The latter can 
be interpreted as the amplitude for the transition of the 
meson with momentum pM into two quarks with mo- 
menta pmx and pm(1 — x). The experimentally derived 
photon-meson transition form factors can be used to test 
different models for the DA. 

The rj and 77' transition form factors have been mea- 
sured in two-photon reactions in several previous exper- 
iments [H-Q. The most precise data for the 77*^'^ at large 
were obtained by the CLEO experiment Q. They 
cover the region from 1.5 to about 20 GeV^. In this 
article we study the 77 and 77' form factors in the range 
from 4 to 40 GeV^. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

We analyze a data sample corresponding to an inte- 
grated luminosity of about 469 fb~^ recorded with the 
BaBAR detector at the PEP-II asymmetric-energy 
storage rings at the SLAG National Accelerator Labo- 
ratory At PEP-II, 9-GeV electrons colhde with 3.1-GeV 
positrons to yield a center-of-mass (cm.) energy near 
10.58 GeV (i.e., the T(4S) resonance peak). About 90% 
of the data used in the present analysis were recorded on- 
resonance and about 10% were recorded about 40 MeV 
below the resonance. 

Gharged-particle tracking is provided by a five-layer 
silicon vertex tracker and a 40-layer drift chamber, op- 
erating in a 1.5-T axial magnetic field. The transverse 
momentum resolution is 0.47% at 1 GeV/c. Energies 
of photons and electrons are measured with a GsI(Tl) 
electromagnetic calorimeter with a resolution of 3% at 
1 GeV. Gharged-particle identification is provided by 
specific ionization (d£'/da;) measurements in the vertex 
tracker and drift chamber and by an internally reflect- 
ing ring-imaging Gherenkov detector. Electron identifi- 
cation also makes use of the shower shape in the calorime- 
ter and the ratio of shower energy to track momentum. 
Muons are identified in the instrumented flux return of 
the solenoid, which consists of iron plates interleaved 
with either resistive plate chambers or streamer tubes. 



Signal e+e~ — )■ e"'"e~77(') and two-photon background 
processes are simulated with the Monte Garlo (MG) event 
generator GGResRc [ll| . It uses the formula for the dif- 
ferential cross section from Ref. Q for pseudoscalar me- 
son production and the Budnev-Ginzburg-Meledin-Serbo 
formalism [l3| for the two-meson final states. Because the 
distribution is peaked near zero, the MG events are 
generated with a restriction on the momentum transfer 
to one of the electrons: > 3 GeV^. This restriction 
corresponds to the limit of detector acceptance for the 
tagged electron. The second electron is required to have 
momentum transfer — g| < 0.6 GeV^. The experimental 
criteria providing these restrictions for data events will 
be described in Sec. IIIII The form factor is fixed to the 
constant value F{Q, 0) in the simulation. 

The GGResRc event generator includes next-to- 
leading-order radiative corrections to the Born cross sec- 
tion calculated according to Ref. [l^. In particular, 
it generates extra soft photons emitted by the initial- 
and final-state electrons. The formulae from Ref. 
are modifled to take into account the hadron contribu- 
tion to the vacuum polarization diagrams. The maxi- 
mum energy of the photon emitted from the initial state 
is restricted by the requirement^ E* < 0.05^/s, where 
■y/s is the e+e^ cm. energy. The generated events are 
subjected to a detailed detector simulation based on 
GEANT4 [Hi and are reconstructed with the software 
chain used for the experimental data. Temporal varia- 
tions in the detector performance and beam background 
conditions are taken into account. 



III. EVENT SELECTION 

The decay modes with two charged particles and two 
photons in the final state, 77' — >■ TT'^Tr~r], 77 — >■ 77 and 
77 — > 7r"'"7r~7r'^, ir^ — >■ 77, are used to reconstruct 77' and 
77 mesons, respectively. For the e^e^ — > e+e~77 process, 
77 — > TT+TT^TT^ is the only decay mode available for anal- 
ysis at BABAR. The trigger efficiency for events with 77 
decays to 2j and to 3tt° is very low. 

Events with at least three charged tracks and two pho- 
tons are selected. Since a significant fraction of signal 
events contains beam-generated spurious track and pho- 
ton candidates, one extra track and any number of extra 
photons are allowed in an event. The tracks correspond- 
ing to the charged pions and electron must have a point 
of closest approach to the nominal interaction point (IP) 
that is within 2.5 cm along the beam axis and less than 
1.5 cm in the transverse plane. The track transverse mo- 
mentum must be greater than 50 MeV/c. The identi- 
fied pion candidates must have polar angles in the range 



^ Throughout this article an asterisk superscript denotes quantities 
in the e+e~ cm. frame. In this frame the positive z-axis is 
defined to coincide with the e~ beam direction. 
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25.8° < 9 < 137.5°, while the track identified as an elec- 
tron must be in the angular range 22.2° < 9 < 137.5° 
(36.7-154.1° in the e+e~ cm. frame). The angular re- 
quirements are needed for good electron and pion identi- 
fication. Electrons and pions are selected using a likeli- 
hood based identification algorithm, which combines the 
measurements of the tracking system, the Cherenkov de- 
tector, and the electromagnetic calorimeter. The elec- 
tron identification efficiency is about 98-99%, with a 
pion-misidentification probability below 10%. The pions 
are identified with about 98% efficiency and a electron- 
misidentification rate of about 7%. To recover electron 
energy loss due to bremsstrahlung, both internal and in 
the detector material before the drift chamber, the energy 
of any calorimeter shower close to the electron direction 
(within 35 and 50 mrad for the polar and azimuthal an- 
gle, respectively) is combined with the measured energy 
of the electron track. The resulting cm. energy of the 
electron candidate must be greater than 1 GeV. 

The photon candidates are required to have laboratory 
energies greater than 50 MeV. For the e+e^ — > e+e"??' 
selection, two photon candidates are combined to form an 
77 candidate. Their invariant mass is required to be in the 
range 0.480-0.600 GcV/c^. To suppress combinatorial 
background from spurious photons, the photon hclicity 
angle is required to satisfy the condition | cos9h\ < 0.9"^. 
The helicity angle 9h is defined in the 77 rest frame as 
the angle between the decay photon momentum and di- 
rection of the boost from the laboratory frame. Each 
candidate is then fit with an 77-mass constraint to im- 
prove the precision of its momentum measurement. An 
77' candidate is formed from a pair of oppositely-charged 
pion candidates and an rj candidate. The 77' invariant 
mass must be in the range 0.920-0.995 GeV/c^. The 77' 
candidate is also then fit with a mass constraint. 

Similar selection criteria are used for e~^e~ — > 6^6^77 
candidates. An 77 candidate is formed from a pair of 
oppositely charged pion candidates and a tt^ candidate, 
which is a combination of two photons with invariant 
mass between 0.115 and 0.150 GeV/c^ and the cosine of 
the photon hclicity angle |cos6';i| < 0.9. The mass of 
the 77 candidate must be in the selection region 0.48-0.62 
GeV/c2. 

Figure [5] shows the | cos 0*^ | distribution for data and 
simulated e'^e~ — > e"'"e~77 events passing the selection 
criteria described above, where 6*1^ is the polar angle of 
the momentum vector of the 677 system in the e"'"e~ cm. 
frame. We require that |cos0*^| be greater than 0.99. 
This condition effectively limits the value of the momen- 
tum transfer to the untagged electron (gf) and guarantees 
compliance with the condition — < 0.6 GeV^ used in 



^ Spurious photons tend to have low energy, and therefore align 
opposite to the rj/ir^ candidate's boost direction, whereas true 
ti/tt^ meson decays into two photons have a flat cos9^ distribu- 
tion. 
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FIG. 2: The |cosSe^| distribution for data events (solid 
histogram). Tlie shaded histogram shows the same distri- 
butions for the e^e~ — e+e~77 simulation. Events with 
I COSSET? I > 0.99 (indicated by the arrow) are retained. 
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FIG. 3: The r distributions for e^e~ — ^ e^e~rj data (solid-line 
histogram) and signal simulation (shaded histogram). The 
arrows indicate the region used to select event candidates 
(-0.025 <r < 0.05). 

the MC simulation. The same condition | cos0*^, | > 0.99 
is used to select the e+e" — ?• e~^e~i]' event candidates. 

The emission of extra photons by the electrons involved 
leads to a difference between the measured and actual 
values of Q^. In the case of initial-state radiation (ISR) 
Qlieas = QlrueiX + ^-y), whcrc = 2E*lsfs- To rcstrict 
the energy of the ISR photon we use the parameter 

J^-E* - \p* „)| 

Ts ' 
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where E*^^,-^ and are the cm. energy and momen- 

tum of the detected ery*-'' system. For ISR this parameter 
coincides with defined above. The r distributions for 



data and simulated e+e 



rj events passing the 



selection criteria described above are shown in Fig. |3l 
For both processes under study, we select events with 
—0.025 < r < 0.05. It should be noted that this condi- 
tion on r ensures compliance with the restriction < 0.1 
used in the simulation. 

For two-photon events with a tagged positron (elec- 
tron), the momentum of the detected erj'^'^ system in the 
e'^e^ cm. frame has a negative (positive) z-component, 
while events resulting from e^e~ annihilation are pro- 
duced symmetrically. To suppress the e^e~ annihilation 
background, event candidates with the wrong sign of the 
momentum z-component are removed. 

The distributions of the invariant masses of rj and rj' 
candidates for data events satisfying the selection criteria 
described above are shown in Fig. |4l For events with 
more than one e^-q^'^ candidate (about 5% of the selected 
events), the candidate with smallest absolute value of the 
parameter r is selected. Only events with 4 < < 40 
GeV^ are included in the spectra of Fig. S) For < 4 
GeV^ the detection efficiency for single-tag two-photon 
77 and 77' events is small (see Sec. IVI[) . In the region 
> 40 GeV^ we do not see evidence of 77 or 77' signal 
over background. About 4350 and 5200 events survive 
the selection described above for 77 and 77', respectively. 



IV. 



FITTING THE Tv+TT-n" AND tt+tt-t? MASS 
SPECTRA 



To determine the number of events containing an 77*^'^ , 
we perform a binned likelihood fit to the spectra shown 
in Fig. 21 with a sum of signal and background distribu- 
tions. The signal distributions are obtained by fitting 
mass spectra for simulated signal events. The obtained 
functions then are modified to take into account a pos- 
sible difference between data and simulation in detector 
response. The signal line shape in simulation is described 
by the following function; 



where 



F{x) = A[G{x) sin^ C + B{x) cos^ C], 



{x - a)' 



G{x) = exp — 



2a2 



B{x) 



(ri/2)'^i 

(x-a)*'2+(r2/2)f2 



if X < a; 
if X > a, 



(3) 



(4) 



(5) 



C, a, cr, Fi, T2, and /32 are resolution function param- 
eters, and A is a normalization factor. The B{x) term 
is added to the Gaussian function to describe the asym- 
metric power-law tails of the detector resolution function. 



The mass spectra for simulated signal events weighted to 
yield the dependencies observed in data and fitted 
curves are shown in Fig. [5l 

When used in data, the parameters cr, Fi, r2 and a 
are modified to account for possible differences between 
data and simulation in resolution (Act) and mass scale 
calibration (Aa): 



'MC 
'MC 



Acr2 if Act < 0; 
Acr2 if Act > 0, 



^ i.MC 
^ j,MC 



(2.35Acr)2 if Act < 0; 
(2.35A(t)2 if Act > 0, 



a = aMC + Aa, 



(6) 



(7) 



(8) 



where the subscript MC indicates the parameter value 
determined from the fit to the simulated mass spectrum. 
The resolution and mass differences. Act and Aa, are 
determined by a fit to data. 

The background distribution is described by a linear 
function. Five parameters are determined in the fit to 
the measured mass spectrum: the number of 77*^'^ events, 
Aa, Act, and two background shape parameters. The fit- 
ted curves are shown in Fig. 2) The numbers of 77 and 77' 
events are found to be 3060 ± 70 and 5010 ± 90, respec- 
tively. The mass shifts are Aa = 0.25 ± 0.09 MeV/c^ for 
the 77 and Aa = -(0.48 ± 0.06) MeV/c^ for the 77'. To 
check possible dependence of the mass shift on , sepa- 
rate fits are performed for two regions: 4 < < 10 
GeV2 and 10 < < 49 QeV^. The Aa values ob- 
tained for these regions agree with each other both for 
ij and rj' . In contrast, the values of Aa are found to 
be strongly dependent on Q^, changing from 0.9 ± 0.3 
MeV/c^ for 4 < < 10 GeV^ to -(1.0 ± 0.6) MeV/c^ 
for 10 < Q2 < 40 GeV^. It should be noted that the 
mass resolution for 77 and 77' is about 4 MeV/c^. The 
data-MC difference, Act 1 MeV/c^, corresponds to a 
small (~ 3%) change in the mass resolution when added 
in quadrature. 

A fitting procedure similar to that described above is 
applied in each of the eleven intervals indicated in 
Table IH The parameters of the mass resolution function 
are taken from the fit to the mass spectrum for simulated 
events in the corresponding interval. The 77 and rj' 
masses are fixed to the values obtained from the fit to 
the spectra of Fig. [D The Act parameter is set to zero. 
Fits with Act = 0.9 MeV/c^ and Act = -1.0 MeV/c^ arc 
also performed. The differences between the results of 
the fits with zero and non-zero Act provide an estimate 
of the systematic uncertainty associated with the data- 
MC simulation difference in the detector mass resolution. 



For the analysis of the e 



'rj process, the 



numbers of events containing an 77 are determined in two 
regions of the parameter r: —0.025 < r < 0.025 (A^i) 
and 0.025 < r < 0.050 {N2). The TVi and N2 values 
are used to determine the numbers of signal events {Ng) 
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FIG. 5: The -k'^-k tt" and 7r"'"7r 77 mass spectra for simulated (a) e^e e^e -q and (b) e^e e^e -q' events, respectively. 
The curves represent the resolution functions described in the text. 



and background events peaking at the 77 mass [Ni,) as de- 
scribed in Sec.|Vl These values arc listed in Table HI The 
Tr+Tr^TT*^ mass spectra and fitted curves for three repre- 
sentative intervals are shown in Fig. [51 The spectra 
shown are obtained for the —0.025 < r < 0.025 regions; 
the 0.025 < r < 0.050 regions contain only 10-13% of 
the signal events and are used mainly to estimate back- 
grounds. 

For the e+e^ — > e+e^?/ process, background is as- 
sumed to be small. There is no need to separate events 
into two r regions. The ■K'^'K~ri mass spectra and fitted 
curves for three representative intervals are shown in 
Fig. [71 The numbers of signal 77' events obtained from 
the fits are listed in Table HIl 



V. PEAKING BACKGROUND ESTIMATION 
AND SUBTRACTION 

Background events containing true 77 or 77' mesons 
might arise from e+e~ annihilation, and two-photon pro- 
cesses with higher multiplicity final states than our sig- 
nal events. The e+e~ annihilation background is studied 
in Sec. IV Al In Sec. IV Bl we use events with an extra 
7r° to estimate the level of the two-photon background 
and study its characteristics. In Sec. IV CI we develop 
a method of background subtraction based on the dif- 
ference in the r distributions for signal and background 
events. This method gives an improvement in accuracy 
compared to the previous one described in Sec. IV Bl and 
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TABLE I: The interval, number of detected e^e~ — ^ e^e~r) signal events (A'^s), number of peaking-background events {N},), 
efficiency correction ((5totai), number of signal events corrected for data-MC difference and resolution effects (A^^o"?''*'^'*), and 
detection efficiency obtained from simulation (e). The first and second errors on Ns and A"o"r°''^'^'* are statistical and systematic, 
respectively. The errors on N}, are statistical and systematic combined in quadrature. 



Q2 interval (GeV^) 


A^. 


Nb 


5totai(%) 


A r unfolded 
' corr 


e(%) 


4-5 


638 ± 31 ± 16 


53 ±27 


-1.4 


634 ± 34 ± 18 


6.3 


5-6 


625 ± 34 ± 19 


89 ±34 


-1.6 


641 ± 38 ± 22 


13.0 


6-8 


622 ± 36 ± 23 


97 ±37 


-1.7 


634 ± 39 ± 25 


14.7 


8-10 


349 ± 26 ± 12 


43 ±23 


-2.0 


359 ± 29 ± 14 


18.7 


10-12 


212 ±20 ±7 


15 ±16 


-2.3 


224 ± 22 ± 8 


22.6 


12-14 


104 ± 14 ± 4 


13± 11 


-2.1 


105 ±17 ±5 


22.9 


14-17 


109 ± 13 ± 3 


0.0 ±9.2 


-2.0 


116 ±15 ±4 


22.2 


17-20 


40.5 ±8.3 ±1.2 


0.7 ±5.6 


-2.3 


41.2 ±9.5 ±1.4 


21.3 


20-25 


32.5 ±7.4 ±0.8 


0.0 ±4.2 


-2.4 


34.4 ± 8.3 ± 0.9 


19.6 


25-30 


13.7 ±5.3 ±0.5 


3.1 ±3.5 


-2.7 


14.2 ±6.0 ±0.6 


18.0 


30-40 


13.0 ±4.8 ±0.3 


0.5 ±3.7 


-2.7 


14.1 ±5.3 ±0.3 


15.7 



TABLE 11: The interval, number of detected r)' signal events (As), number of peaking-background events (A(,), efhciency 
correction (5totai), number of signal events corrected for data-MC difference and resolution effects (A^q"'°''*°'^), and detection 
efficiency obtained from simulation (e). The first and second errors on A^s and A"orr°''*'^'* are statistical and systematic, 
respectively. 



interval (GeV^) 


A^. 


Nt, 


5total(%) 


Tvrunfoldcd 
*corr 


e(%) 


4-5 


950 ± 32 ± 5 


0.0 ±0.0 


-0.4 


936 ± 34 ± 6 


5.7 


5-6 


1013 ± 33 ± 6 


0.0 ±0.0 


-0.6 


1015 ± 36 ± 7 


12.5 


6-8 


1185 ± 36 ± 5 


0.0 ±0.0 


-0.7 


1207 ± 38 ± 6 


14.3 


8-10 


710 ± 28 ± 3 


0.0 ±0.0 


-1.0 


716 ± 30 ± 4 


19.9 


10-12 


454 ± 22 ± 4 


0.0 ±0.0 


-1.2 


467 ± 25 ± 4 


26.4 


12-14 


243 ± 16 ± 1 


0.0 ±0.0 


-1.0 


250 ± 19 ± 1 


28.1 


14-17 


207 ± 15 ± 2 


0.0 ±0.0 


-0.8 


214 ±17 ±2 


28.1 


17-20 


108 ± 10 ± 1 


0.0 ±0.0 


-0.8 


112 ± 12 ± 1 


26.8 


20-25 


80.0 ±9.0 ±0.1 


0.0 ±0.0 


-1.0 


82.5 ± 9.9 ± 0.2 


26.3 


25-30 


30.2 ±5.9 ±0.2 


1.0 ± 1.0 


-1.3 


31.7 ±6.7 ±0.2 


25.6 


30-40 


17.2 ±5.4 ±0.1 


2.0 ± 1.4 


-1.4 


18.1 ± 5.8 ± 0.1 


22.5 



has a lower sensitivity to the model used for background 
simulation. 



A. e e annihilation background 

The baekground from e^e~ annihilation ean be esti- 
mated using events with the wrong sign of the e^rj'^'^ 
momentum z-component. The numbers of background 
events from e^e~ annihilation in the wrong- and right- 
sign data samples are expected to be approximately 
the same, but their distributions are quite different. 
The distribution expected for right-sign background 
events coincides with the distribution for wrong-sign 
events, where is the squared difference between the 
four- momenta of the detected positron (electron) and the 
initial electron (positron). 

In the region from 4 to 40 GeV^ we observe 3 
wrong-sign events in the rj' data sample, all peaking at 
the 77' mass, and 9 events in the rj data sample, 5 of 
which are in the 0.530-0.565 GeV/c^ mass window. The 
contribution from non-?7 events to this mass window is 



estimated to be 0.3 events. A possible source of these 
events is the e^e" — > X7 process, where X is a hadronic 
system containing an jy or jy' meson, for example, tt'^tt^t]' , 
with the photon emitted along the beam axis. 

The distribution for the wrong-sign events is used 
to estimate the distribution for e^e^ annihilation 
background in the right-sign data sample. The fraction of 
e+e~ annihilation events in the 77^'' data sample is about 
lO^'^. However, such events are the main contribution 
to the peaking background in high bins and cannot 
be neglected. For the e+e~ — >■ e^e~T]' process, for which 
we do not observe a significant two-photon background 
(see Sec lVBp . the 3 background events from e^e~ anni- 
hilation are subtracted from the two highest intervals 
(see Table HD. 

For the e+e^ — > e^e^rj process, the e'^e" annihilation 
events are effectively subtracted with the procedure de- 
veloped for subtraction of two-photon background (see 
Sec IV C[) . The procedure exploits the difference between 
the r distributions for signal and background events. The 
r distribution for the e'^e~ annihilation events (3 of 5 
events have r > 0.025) is close to that for two-photon 
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FIG. 6: The tt'^-k'tt" mass spectra for data events with 
—0.025 < r < 0.025 for three representative intervals. The 
sohd curves are the fit results. The dashed curves represent 
non-peaking background. 



background. 

In future high statistics measurements of the meson- 
photon form factors at Super B factories e+e~ annihi- 
lation will be the dominant background in the high 
region (Q^ > 50 GeV^). 



B. Two-photon background 



Other possible sources of peaking background are the 
two-photon processes e+e~ — e+e^T^^'-'Tr*'. For the rj 
selection the additional background comes from the two- 
photon production of rj' mesons followed by the decay 
chain r/ — > ir'^n'^ri, rj — > tt+tt^tt*'. The distribution 
of events from the latter background source is calculated 
from the distribution of the selected 77' events. The ra- 



tio of the detection efSeicncies for the two r]' decay modes 
is obtained from MC simulation. The total number of 
77' — > 7r°7r°77 events in the 77 data sample is estimated to 
be 17 ± 2. The events are concentrated almost entirely 
in the three lowest bins. 

To estimate background contributions from the 
e^e~ — >■ e~ ■q^''^ processes, we select events with two 
extra photons that each have an energy greater than 70 
MeV. The distributions of the invariant mass of these 
extra photons for 77 and 77' events are shown in Fig. [S] 
The invariant masses of the 77 and 77' candidates are re- 
quired to be in the mass windows 0.530-0.565 GeV/c^ 
and 0.945-0.970 GeV/c^, respectively. The spectra are 
fit by a sum of the 7r° line shape obtained from simulated 
e+e^ — > e~ r]^'^ TT^ events and a quadratic polynomial. 
The fitted numbers of events with an extra 7r° are 90± 20 
and 13 ± 14 for the 77 and 77' selections, respectively. It is 
expected that 8 events with an extra tt" in the 77 sample 
arise from two-photon 7/ production. 

The distribution of the 777r'^ invariant mass for events 
with an extra is shown in Fig. [S] The two-photon in- 
variant mass of the tt*^ candidate is required to be in the 
0.115-0.150 GeV/c2 range. The sidebands, 0.065-0.100 
and 0.170-0.205 GeV/c^, are used to subtract contami- 
nation from non-777r'^ events. It is known from two-photon 
measurements in the no-tag mode [l^ that the rjn^ final 
state is produced mainly via ao(980) and a2(1320) inter- 
mediate resonances. Evidence for these two intermediate 
resonances is seen in the mass spectrum of Fig. [5) Our 
spectrum differs significantly from the spectrum for the 
no-tag mode [lB|, which is dominated by a2(1320) pro- 
duction. In the no-tag mode the a2(1320) meson is pro- 
duced predominantly in a helicity-2 state, and thus with 
an angular distribution proportional to sin* 6'^ , where 
is the angle between the tt^ direction and the 77 collision 
axis in the 77 cm. frame. Our selection criteria favor 
events with values of 6,^ near zero and hence suppress 
hehcity-2 states. 

From MC simulation we estimate that the ratio of 
the number of e~^e~ — >■ e'^ e~ r/^''^ tt^ events with a de- 
tected TT^ to the number selected with standard criteria 
is about 2.5. For the e+e" — > e^e^Vj process the esti- 
mated two-photon background does not exceed 1.6% of 
the total number of selected 77' events at 90% confidence 
level. This background level is treated as a measure of 
the systematic uncertainty due to possible two-photon 
background for the e+e~ err]' process. 

For the e+e^ — > e^e^rj process the two-photon back- 
ground is about 10% of the total number of selected 77 
events. It should be noted that in the CLEO publica- 
tion on measurements of the meson-photon transition 
form factors the background from the two-photon pro- 
duction of the 7;7r'' final state was not considered. 

A similar technique is used to estimate background 
from the process e~^e~ — e"'"e~0, 4> ^ We do not 
see any (j) meson signal in the 7^7 mass spectrum and 
estimate that this background does not exceed 10% of 
the 7777*^ background. The 7^7 events have the r distri- 
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FIG. 8: The two-photon invariant mass spectra for (a) 77 and (b) 77' events with two extra photons. The sohd histograms 
represent the fit results. The dashed curves are the fitted distributions for events without an extra tv". 



bution similar to that for 7777° events, and are effectively 
subtracted by the procedure described in the next sec- 
tion. The background contributions from the processes 
e"^e~ — ^ e~^e~4>, </> — >■ is negligible due to the small 
(j) — ^ rj' gamma branching fraction. The background from 
e+e~ e^e~J/^j, J/il) ^ 7/(''7 is estimated using the 
distribution of e"*'e^ — )■ e'^e~J/%l} events measured 
in Ref. and efficiencies from MC simulations, and is 
found to be negligible. 



C. Background subtraction from the 77 data sample 

To subtract background from the 77 data sample the 
difference between the r distributions for signal and back- 
ground events is used. The parameter r is proportional 
to the difference between the energy and the momentum 



of particles recoiling against the 677^'^ system and, there- 
fore, is close to zero for signal and has nonzero positive 
value for background events. To obtain the r distribu- 
tion, data events are divided into 15 r intervals. For each 
interval, the fit to the 7r"'"7r~7r'^ (7r"'"7r~77) spectra is per- 
formed and the number of events containing an rj'^'^ is 
determined. The r distributions for events in the 77 and 
rj' data samples are shown in Fig. [TOl 

For 77' events, for which the background is small, the 
data distribution is compared with the simulated signal 
distribution normalized to the number of data events. 
The distributions are in reasonable agreement. The ratio 
Rs of the number of events with r > 0.025 to the number 
with r < 0.025 is found to be 0.103 ± 0.006 in data and 
0.116 ± 0.002 in simulation; the 13% difference is taken 
as a systematic uncertainty on the Rs value for 77' events 
determined from simulation. Since the simulated r dis- 
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based on its rjir'^ mass dependence. The maximum devi- 
ation from the value averaged over the ijir^ spectrum of 
about 25% is found when we exclude events with mass 
near the jyTr" threshold. This deviation is taken as an 
estimate of the systematic uncertainty on Rf,. The r 
distribution for background events from two-photon rj' 
production (i?b is about 10) differs significantly from the 
distribution for tjtt'^ events. Therefore we first subtract 
the calculated r]' contribution from A^i and N2 in each 
interval, and then calculate Ng assuming that the re- 
maining background comes from the e~^e~ — > e'^e~r]TT^ 
process. The obtained numbers of signal and background 
events are listed in Table H] The background includes 
both the e'^e^ — > e+e^?/7r° and e+e~ — >■ e^e~?/ contri- 
butions. The systematic errors quoted for Ng arc mainly 
due to the uncertainties on Rg and Rb- 



FIG. 9: The distribution of the rjiv invariant mass for rj events 
with an extra tt". The background from non-r;7r° events is 
subtracted. 

tributions for 77' and rj events are very close, the same 
systematic error can be applied to Rg value for 77 events. 

For 77 events, the data r distribution is fit with the sum 
of the simulated distributions for signal and background 
e+e~ — > e^e~r]Tr^ and e+e~ — > e^e^ri' — > e^e~Tr^Tr^ri 
events. The fitted number of background events is 280 ± 
40, in reasonable agreement with the estimate given in 
the previous subsection based on the number of events 
with a detected extra tt". 

To subtract the background in each interval the 
following procedure is used. In Sec. llVl we described how 
the number of events containing an rj is determined for 
two regions of the parameter r: —0.025 < r < 0.025 (A^i) 
and 0.025 < r < 0.050 {N2). The numbers of signal and 
background events are then calculated as follows: 

_ {1 + Rs){NiRb ~ N2) 



«.^<i±Mi!j_^, (10) 

where Rg {Rb) is the N2/N1 ratio obtained from sig- 
nal (background) MC simulation. The expressions in 
Eqs. © and PH]) are equivalent to a two-r-bin fit of 
data to signal and background MC predictions; fits us- 
ing a higher number of bins are not useful due to lack of 
statistics. 

The parameter Rg is found to vary from 0.15 to 0.10 
with increasing Q^. The systematic uncertainty on Rg 
(13%) was estimated above. To calculate Rb for the 
e~^e~ — >■ e^e^TyTr" process, the simulated background 
events are reweighted to reproduce the jyTr" mass spec- 
trum observed in data (Fig.lH]). The Rb value varies from 
2.0 to 1.5. The systematic uncertainty on Rb is estimated 



VI. DETECTION EFFICIENCY 

The detection efficiency is determined from MC simu- 
lation as the ratio of the true distributions computed 
after and before applying the selection criteria. The 
dependencies of the detection efficiencies for both pro- 
cesses under study are shown in Fig. [11] The detector ac- 
ceptance limits the detection efficiency at small Q^. The 
cross sections are measured in the regions > 4 GeV^, 
where the detection efficiencies are greater than 5%. The 
asymmetry of the e'*"e~ collisions at PEP-II leads to dif- 
ferent efficiencies for events with electron and positron 
tags. The range from 4 to 6 GeV^ is measured only 
with the positron tag. 

We study possible sources of systematic uncertainty 
due to differences between data and MC simulation in 
detector response. The MC simulation predicts about a 
2.5% loss of signal events, weakly dependent on Q"^, due 
to the offline trigger, i.e. program filters, which provide 
background suppression before the full event reconstruc- 
tion. Events of the process under study satisfying our se- 
lection criteria pass a filter selecting events with at least 
three tracks in the drift chamber originating from the in- 
teraction region. The filter inefficiency is measured from 
data using a small fraction of selected events that does 
not pass the background filters. Combining events from 
the 77 and rj' samples, we determine the ratio of the in- 
efficiencies in data and MC simulation to be 1.15 ± 0.20. 
The error of the ratio is used to estimate the systematic 
uncertainty for the filter inefficiency: 0.2 x 2.5 = 0.5%. 
The trigger inefficiency obtained using MC simulation is 
about 1% in the first interval (4-5 GeV^) and falls 
to zero at > 14 GeV^. The limited statistics do not 
allow us to measure this inefficiency in data. Therefore, 
the level of the inefficiency observed in the MC simula- 
tion is taken as an estimate of the systematic uncertainty 
due to the trigger inefficiency. 

The systematic uncertainty due to a possible difference 
between data and simulation in the charged-particle track 
reconstruction for pions is estimated to be about 0.35% 
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FIG. 10: (a) The r distribution for data events containing an ry (points with error bars). The dashed histogram shows the 
fit results. The shaded histogram is the fitted background contribution from the processes e^e" e^e^-q-K^ and e'" — >■ 
e^e~r)' — >■ e~ 'K^iv'^r). (b) The r distribution for data events containing an r)' (points with error bars). The solid histogram is 

e~ri 



the simulated distribution for events from the signal e e 



process normalized to the number of data events. 
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FIG. 11: The detection efficiencies for (a) e^e~ — >■ e^e~r] with rj — >■ ■K'^ir~-n^ and tt" 27 and (b) e^e~ e^e~r] with 
77 — 7- ■K'^Tv~rj and 7; — >■ 27 as functions of the momentum transfer squared for events with a tagged electron (squares), a tagged 
positron (triangles), and their sum (circles). In the region < 6 GeV^, where the electron-tag efficiency is close to zero, the 
sum and the positron-tag efficiencies coincide. 



per track, so the total uncertainty is 0.7%. For electron 
tracks, this uncertainty is about 0.1%. 

The data-MC simulation difference in the pion identi- 
fication efficiency is estimated using the identification ef- 
ficiencies measured for pions in the Z?*+ — D'^tt^, — > 
Ti^K~ decay. The ratio of the data and MC identifica- 
tion efliciencies is determined as a function of the pion 
momentum and polar angle. These functions for positive 
and negative pions are then convolved with the pion en- 
ergy and angular distributions for simulated signal events 



in each interval. The resulting efficiency correction 
((Jtt) for pion identification varies from -1% to 0.5% in 
the range from 4 to 40 GeV^. The systematic un- 
certainty in the correction does not exceed 0.5%. The 
data-MC simulation difference in electron identification 
is estimated using the identification efficiencies measured 
for electrons in radiative Bhabha events. The found effi- 
ciency correction {5^) does not exceed 1%. Its systematic 
uncertainty is estimated to be 0.5%. 

The reconstruction efficiency is studied using events 
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FIG. 12: The correction to the MC-estimated 7r° reconstruc- 
tion eSiciency 5^o as a function of for the e'^e~ — > e'^e~ri 
process. 
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from the ISR process e+e" — 7W, u 7r+7r~7r°. These 
events can be reconstructed and selected without using 
information related to the tt*^. The tt*^ reconstruction ef- 
ficiency is computed as the ratio of the number of events 
with an identified tt*^ to the total number of reconstructed 
e+e^ — > 7Cl> events. The data-MC simulation relative 
difference in the 7r° efficiency depends on the tt^ momen- 
tum and varies from (0.7± 1.2)% at momenta below 0.25 
GeV/c to (-4.2 ± 1.3)% at 4 GeV/c ^. The eflnciency 
correction averaged over the ir^ spectrum is shown in 
Fig. HHas a function of Q^. The systematic uncertainty 
associated with this correction is estimated to be 1%. 
For 77 — >• 77 decays the efficiency correction is expected 
to be smaller. The maximum value of the 7r° efficiency 
correction (2%) is conservatively taken as an estimate of 
systematic uncertainty due to a possible data-MC simu- 
lation difference in the 77 — >■ 77 reconstruction. 

To estimate the effect of the requirement —0.025 < r < 
0.05, 77' events with 0.05 < r < 0.075 are studied. We 
calculate the double ratio minus unity 

Act _ (iVncw/A^)data ^ 
(J (iV„ew/iV)MC ' ^ ' 

where Nnew and N are the numbers of signal events with 
the new and standard selection criteria. The ratio is sen- 
sitive to the relative change in the measured cross sec- 
tion due to the changes in the selection criteria. Wc 
do not observe any significant dependence of Aa ja. 
The average over is found to be consistent with zero 
(-0.003 ± 0.004). We conclude that the simulation re- 
produces the shape of the r distribution. 

We also study the effect of the |cos0*^,,)| > 0.99 re- 
striction by changing the value to 0.95. The correspond- 
ing change of the measured cross section does not depend 
on . The average change in cross section integrating 
over is (2.0 ± 0.4)%. We consider this data-MG sim- 
ulation difference (2%) as a measure of the systematic 
uncertainty due to the cos9*^^f^,^ criterion. 

The angular and energy distributions of detected par- 
ticles are very different for events with electron and 



FIG. 13: The ratio of the cross sections for the e e~ — ^ e e~?7 
process measured with electron and positron tags as a func- 
tion of Q^. 

positron tags. As a cross-check of our study of the ef- 
ficiency corrections, we have performed comparison of 
dependencies of the cross sections obtained with only 
electron and only positron tags. For > 8 GeV^, where 
both positron and electron data are available, the ratio of 
the cross sections have been found to be consistent with 
unity, for both rj and rj' events. The dependence of 
the ratio for ry events is shown in Fig. 1131 Due to limited 
statistics data of the three highest bins are combined. 

The main sources of systematic uncertainty associ- 
ated with the detection efficiency are summarized in Ta- 
ble mil for both processes under study. The values of 
the detection efficiency and the total efficiency correc- 
tion (5totai = + 6e + ^7r" (thc tcmi S^^o IS Only applicable 
to the ?7 mode) for different intervals are listed in 
Tables U and [TTl The data distribution is corrected as 
follows: 

Nr"- ^Nj{l + Stot.L,), (12) 

where Ni is the number of signal events in the ith 
interval. 

VII. CROSS SECTION AND FORM FACTOR 

The Born differential cross section for e+e^ — > 
e^e~T]^'^ is 

dg2 eRLB ^ ' 

where (dA^/dQ^)™/"'^'"^ is thc mass spectrum corrected 
for data-MC simulation differences and unfolded for de- 
tector resolution effects as explained below, L is the to- 
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TABLE III: The main sources of systematic uncertainty asso- 



ciated with the detection efficiency, 


and the 


total efficiency 


systematic uncertainty for e^e~ — 




and f 




e^e~ri' events. 










n (%) 




v' (%) 


Track reconstruction 




0.8 




TT* identification 




0.5 




identification 




0.5 




|cosS*^(,)| > 0.99 criterion 
Trigger, filters 




2.0 




0.7 




1.3 


77, -K^ 27 reconstruction 


1.0 




2.0 


Total 


2.6 




3.3 



tal integrated luminosity, e is the Q^-dependent detec- 
tion efficiency, and i? is a radiative correction factor ac- 
counting for distortion of the spectrum due to vac- 
uum polarization effects and the emission of soft photons 
from the initial-state particles. The factor B is the prod- 
uct of the branching fractions, B{ri — > 7r+7r^7r")S(7r" 
77) = 0.2246 ± 0.0028 or B{tj' Tr+Tr-r])B{r] 77) = 
0.1753 ±0.0056 0. 

The radiative correction factor R is determined using 
simulation at the generator level, i.e., without detector 
simulation. The spectrum is generated using only 
the pure Born amplitude for the e~^e~ — > e^e~r]^'^ pro- 
cess, and then using a model with radiative corrections 
included. The radiative correction factor, evaluated as 
the ratio of the second spectrum to the first, varies from 
0.994 at Q2 = 4 GeV^ to 1.002 at = 40 GeV^. The 
accuracy of the radiative correction calculation is esti- 
mated to be 1% [ll|. It should be noted that the value 
of R depends on the requirement on the extra photon 
energy. The dependence obtained corresponds to the 
condition r — 2E*/y/s < 0.1 imposed in the simulation. 

The corrected and unfolded distribution 
(diV/dQ2)unfoidcd ig obtained from the measured 

distribution by dividing by the efficiency correction 
factor (see Eq. p^ ) and unfolding for the effect of finite 
resolution. Using MC simulation, a migration matrix 
H is obtained, which represents the probability that 
an event with true in interval j is reconstructed in 
interval i: 




In the case of extra photon emission, Q^^.^^^ is calculated 
as —{p — p' — k)^, where k is the photon four-momentum; 
e and R in Eq. (fT3)) are functions of Qtmc- the cho- 
sen interval width significantly exceeds the resolution 
for all Q^, non-zero elements of the migration matrix lie 
on and near the diagonal. The values of the diagonal ele- 
ments are in the range 0.9-0.95. The true distribution 
is obtained by applying the inverse of the migration ma- 
trix to the measured distribution. The procedure docs 
not change the shape of the distribution significantly. 



but increases the errors (by about 10%) and their corre- 
lations. The number of events (A^^"„';.f°i'^°^) as a function 
of is reported in Tables U and HIl 

The value of the differential cross section as a func- 
tion of is listed in Tables |TV] and El The quoted 
errors are statistical and systematic. The latter includes 
only (5^-dependent errors: the systematic uncertainty in 
the number of signal events and the statistical errors on 
the efficiency correction and MC simulation. The Q^- 
independent systematic error on the e~^e~ — >■ e^e~rj cross 
section is 3.5%; this results from the uncertainties on the 
detection efficiency, both systematic (2.6%) and model- 
dependent (1.5%), the uncertainty in the calculation of 
the radiative correction factor (1%), and the errors on 
the integrated luminosity (1%) and the 77 decay branch- 
ing fraction (1.2%) [l^. The Q^-independent systematic 
error on the e+e~ — > e^e~rj' cross section is 5.3%. It 
includes the systematic and model uncertainties on the 
detection efficiency (3.3% and 1.5%, respectively), the 
uncertainties on the background subtraction (1.6%) and 
the radiative correction factor (1%), and the errors on 
the integrated luminosity (1%) and the rj' decay branch- 
ing fraction (3.2%) 0. 

The model dependence of the detection efficiency arises 
from the unknown cross-section dependence on the mo- 
mentum transfer to the untagged electron. The MC sim- 
ulation is performed, and the detection efficiency is de- 
termined, with the restriction that the momentum trans- 
fer to the untagged electron be greater than —0.6 GcV^, 
so that the cross section is measured for the restricted 
range jgH < 0.6 GeV^. The actual (7^ threshold is de- 
termined by the requirement on cos and is equal to 

0.38 GeV^. The MC simulation is performed with a 
independent form factor, which corresponds to the QCD- 
inspired model F{qf,q'^) oc l/{qf + g|) « 1/qf [l7|. The 
event loss due to the jgH < 0.38 GeV^ restriction is about 
2.5%. The use of the form factor predicted by the vec- 
tor dominance model F(g|) oc 1/(1 — g|/m^), where uip 
is p meson mass, leads to a decreased event loss of only 
1%. The difference between these efficiencies is consid- 
ered to be an estimate of the model uncertainty due to 
the unknown g| dependence. 

Because of the strong nonlinear dependence of the 
cross section on Q^, the effective value of correspond- 
ing to the measured cross section differs from the center 
of the interval. We parametrize the measured cross 
section with a smooth function and calculate for each 
interval solving the equation 

da/d(g2)(Q^) = d<j/d{Q\^,,,^,, 

where d(T/d((5^)avorago is the differential cross section av- 
eraged over the interval. The values of are listed 
in Table IIVI and Table |Vl The measured differential 
cross sections for both processes under study are shown 
in Fig. [T31 together with the data reported by the CLEO 
Collaboration Q for > 3.5 GeV^. We average the 
CLEO results obtained in different 77*^'' decay modes as- 
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TABLE IV: The interval, the weighted average Q value 
for the interval (Q^), the e^e~ — e^e~ri cross section 
(dcr/dQ'^(Q2)), and the product of the 77* — > r) transition 
form factor F{Q^) and Q^. The statistical and systematic 
errors are quoted separately for the cross sections, and are 
combined for the form factors. In the table we quote the 
Q'^-dependent systematic errors. The Q^-independent error 
is 3.5% for the cross section and 2.9% for the form factor. 



interval 
(GeV") 


(GeV") 


da/dQ2(Q2) 
(fb/GeV^) 


Q'\F{Q^)\ 

(MeV) 


4-5 


4.47 


95.6 ±5.1 ±3.1 


143.4 ±4.4 


5-6 


5.47 


46.6 ± 2.7 ± 1.7 


142.7 ±4.9 


6-8 


6.89 


20.4 ±1.2 ±0.8 


142.6 ± 5.2 


8-10 


8.92 


9.06 ±0.72 ±0.35 


151.2 ±6.7 


10-12 


10.96 


4.67 ±0.47 ±0.18 


158.5 ±8.5 


12-14 


12.92 


2.16 ±0.34 ±0.10 


146.5 ± 12.1 


14-17 


15.38 


1.65 ±0.22 ±0.06 


178.9 ± 12.1 


17-20 


18.34 


0.61 ±0.14 ±0.02 


151.6 ±17.8 


20-25 


22.33 


0.33 ±0.08 ±0.01 


166.0 ± 20.2 


25-30 


27.23 


0.15 ±0.06 ±0.01 


166.7 ±36.6 


30-40 


34.38 


0.085 ± 0.032 ±0.003 


205.9 ±39.0 



TABLE V: The interval, the weighted average value 
for the interval (Q^), the e'^e~ — > e~ r] cross section 
(d(T/dQ^((3^)), and the product of the 77* — ^ 77' transition 
form factor F{Q'^) and Q^. The statistical and systematic 
errors are quoted separately for the cross sections, and are 
combined for the form factors. In the table we quote the 
Q'^-dependent systematic errors. The Q^-independent error 
is 5.3% for the cross section and 3.5% for the form factor. 



interval 
(GeV') 


(GeV") 


da/dQ2(Q2) 
(fb/GeV^) 


Q^\Fm\ 

(MeV) 


4-5 


4.48 


202 ± 7 ± 3 


216.2 ±4.3 


5-6 


5.46 


99.6 ±3.6 ±1.4 


214.3 ±4.1 


6-8 


6.90 


51.7 ± 1.6 ±0.5 


233.3 ±3.9 


8-10 


8.92 


22.1 ±0.9 ±0.2 


241.6 ±5.2 


10-12 


10.95 


10.8 ±0.6 ±0.1 


245.5 ±6.7 


12-14 


12.90 


5.45 ±0.41 ±0.06 


236.7 ±8.9 


14-17 


15.33 


3.10 ±0.24 ±0.04 


248.5 ±9.9 


17-20 


18.33 


1.70 ±0.18 ±0.02 


258.7 ± 13.7 


20-25 


22.36 


0.77 ±0.09 ±0.01 


257.0 ± 15.4 


25-30 


27.20 


0.30 ±0.06 ±0.01 


240.0 ± 25.7 


30-40 


34.32 


0.098 ±0.031 ±0.002 


224.1 ±35.9 



suming that systematic errors for different modes are not 
correlated. 

To extract the transition form factor, the measured 
and calculated cross sections are compared. The simu- 
lation uses a constant form factor F^^. Therefore, the 
measured form factor is determined from 



\Fm\' 



data 

(da/dQ2) 



MC 



MC- 



(15) 



The calculated cross section {da / dQ^)Mc has a model- 
dependent uncertainty due to the unknown dependence 
on the momentum transfer to the untagged electron. 



The difference between the cross section values calcu- 
lated with the two form-factor models described above 
is 4.6% for both r/ and 77'. This difference is considered 
to be an estimate of the model uncertainty due to the 
unknown q| dependence. The values of the form factors 
obtained, represented in the form Q'^\F{Q^)\, are listed 
in Tables IIVI and |V] and shown in Fig. [151 For the form 
factor we quote the combined error, obtained by adding 
the statistical and Q^-dependent systematic uncertain- 
ties in quadrature. The Q^-independcnt systematic error 
is 2.9% for the 77 and 3.5% for the 77' form factor. 



VIII. DISCUSSION AND SUMMARY 

The comparison of our results on the form factors with 
the most precise previous measurements @ is shown in 
Fig. [ini For the 77' form factor our results are in good 
agreement with those reported by the CLEO Collabora- 
tion Q. For the 77 form factor the agreement is worse. In 
particular, the CLEO point at w 7 GeV^ lies higher 
than our measurements by about 3 standard deviations. 

The data for the e+e~ — > ?7'^''7 reactions are used to 
determine the transition form factors in the time-like re- 
gion = s > 0. Since the time- and space-like form 
factors are expected to be similar at high Q^, in Fig. \W\ 
we show the results of the high-Q^ time-like measure- 
ments together with the space-like data. The form fac- 
tors at = 14.2 GeV^ are obtained from the values of 
the e'^e~ — >■ 77*^'' 7 cross sections measured by CLEO [l^ 



near the peak of the -0(3770) resonance. We calculate the 
form factor using the formulas from ReL [l^ under the 
assumption that the contributions of the -0(3770) — >■ 7;'^')7 
decays to the e+e^ — > 77(''7 cross sections are negligible. 
It is seen that the measured time- and space-like form fac- 
tors at 14 GeV^ are in agreement both for 77 and 
for ?]' . The BaBAR measurements of the e+e^ rj^'^j 
cross sections [l^ allow us to extend the region for 
the 77 and 77' form factor measurements up to 112 GeV^. 



In most models for the meson distribution amplitude 
4>p{x) used for calculation of photon-meson transition 
form factors, the DA end-point behavior is determined by 
the factor x{l — x). The form factors calculated with such 
conventional DAs are almost flat for values greater 
than 15 GeV^ (see, for example, the recent works [20I - 
[2^ devoted to the 77* — )- 7r° form factor). Some of these 
models [l^l have difficulties in reproducing the depen- 
dence of the 77* — > 7r° form factor measured by BaBAR [i| 
in the range from 4 to 40 GeV^. Alternatively, mod- 
els with a flat DA or a DA that is finite at the end points 
have been suggested [23l - l25j . which give a logarithmic 
rise of the product Q^F{Q'^) with and describe the 
BaBAR data reasonably well. 

The dependencies of the products Q^F{Q^) for 77 
and 77' are fit with the function 



'F{Q^)^bi + ai lnQ2(GeV2). 
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FIG. 14: The differential cross sections for (a) e^e e^e rj and (b) e e e^e r]' from the present analysis compared to 
those from the CLEO experiment [^. The asterisk near the label "CLEO" in this and next figures indicates that the original 
CLEO results obtained in different r;''' decay modes were averaged assuming that systematic errors for different modes are not 
correlated. In the present analysis the cross sections are measured with the restriction \q2\ < 0.6 GeV'^. In the CLEO analysis 
the cross sections have been obtained using the vector dominance model for the gl dependence in simulation. 
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FIG. 15: The transition form factors multiplied by for (a) 77* — ^ -q and (b) 77* — >■ rj' 



The results of the fit are shown in Fig. [161 For both jy and 
77' the quality of the fit is acceptable: /v is equal to 
6.8/10 for 77 and 15.9/10 for 77', where v is the number of 
degrees of freedom. The observed rise of the form factors 
{ai « 0.20 ± 0.05 GeV) is about three times weaker than 
the corresponding rise of the tt'^ form factor . 

The dashed horizontal lines in Fig. [TCI show the results 



of fits assuming Q^F{Q^) to be constant for 14 < < 
112 GeV^. The average values of Q^F{Q^) in this range 
are 0.175 ± 0.008 GcV for ij and 0.251 ± 0.006 GeV for 
rj' . The x^/v for the fits are 5.6/5 for the rj and 1.3/5 
for the Tj' . The preferred description for the i] form fac- 
tor is the logarithmic function of Eq. p^ . corresponding 
to the models with a finite DA at the end points. The 
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FIG. 16: The transition form factors multiplied by for (a) 77* — !> 77 and (b) 77* — !• 77'. The solid line shows the result of the 
fit to BABAR data by the function given by Eq. (|16|) . The dashed lines indicate the average form factor values over the data 
points with > 14 GeV^ 



77' form factor is better described by the model with a 
conventional DA, yielding a flat Q^F{Q'^) for > 15 
GeV2. 

To compare the measured values of the ?/ and ?/' form 
factors with theoretical predictions and data for the 7r° 
form factor we use the description of 77-?/ mixing in the 
quark flavor basis [2^ : 

I") = ■^{\nu) + \dd)), \s) = \ss), 

\ri) = cos |7i) — sin |s), 
It;') = sin(?!)|7i) +cos0|s), (17) 

where (p is the mixing angle. The 77 and 77' transition form 
factors are related to the form factors for the \n) and |s) 
states: 

= cos F„ — sin Fs , F,,/ = sin F„ + cos , (18) 
which have asymptotic limits for — > 00 [27j given by 

Q'F,{Q^) = lf,, Q^F„iQ')=^-^U, (19) 

where /„ and fs are the decay constants for the |n) 
and \s) states, respectively. For the tt^ form factor, the 
corresponding asymptotic value is v^/tt- The pion de- 
cay constant is determined from leptonic tt decays to be 
130.4 ± 0.2 MeV For the \n) and \s) states, we 

use the "theoretical" values from Ref. [1^: /„ = and 
= V2/i-/2 « 1.36/^ {/k/L = 1.193 ±0.006 0), 
which agree to within 10% with the "phenomenological" 
values [26[ extracted from the analysis of experimental 



data, for example, for the two-photon ?/ and 77' decays. 
The currently accepted value of the mixing angle is 
about 41° [2^. Under the assumption that the \n) and 
TT*^ distribution amplitudes arc similar to each other, the 
only difference between the \n) and tt*' form factors is 
a factor of 3/5 that arises from the quark charges. In 
Fig. [17] the form factor for the |7i) -state multiplied by 
3(5^/5 is compared with the measured 7*7 — > tt*^ form 
factor [l| and the results of the QCD calculations per- 
formed by A. P. Bakulev, S. V. Mikhailov and N. G. Ste- 
fanis [23] for the asymptotic DA [3^, the Chernyak- 
Zhitnitsky 71° DA [HI, and the 7r° DA derived from QCD 
sum rules with non-local condensates [32j . The horizon- 
tal dashed line indicates the asymptotic limit for the tt*^ 
form factor. 

The dependencies of the measured |7t,) and 7r° form 
factors are significantly different. This indicates that the 
distribution amplitudes for the \n) and tt*' are signifi- 
cantly different as well. The data for the \n) form factor 
are well described by the model with DA from Ref. [s^] , 
while the data for the 7r° form factor is reproduced by 
the models with a significantly wider DA [2ll [2^ or a 
flat DA [il-il. 

The form factor for the \s) state is shown in Fig. [151 
The dotted curve shows the QCD prediction [2§| for the 
asymptotic DA (soj . defined by multiplying the tt^ curve 
in Fig. [TTjby a factor of {y/2/2>)fs/f-n- The data lie sys- 
tematically below this prediction. This may indicate, 
in particular, that the distribution amplitude for the \s) 
state is narrower than the asymptotic DA. However, due 
to the strong sensitivity of the result for the the \s) state 
to mixing parameters, other interpretations are possible. 
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For example, an admixture of the two-gluon component 
in the rj' meson (33l - l36j can lead to a significant shift of 
the values of the |s) form factor. 

In summary, we have studied the e^e^ — > e^e~ri and 
e+e~ — > e~^e~ri' reactions and measured the differential 
cross sections (dcr/dQ^) and the 7*7 — > transition 
form factors F{Q^) in the momentum transfer range from 
4 to 40 GeV^. In general, our results are in reasonable 
agreement with the previous CLEO measurements 0]. 
We significantly improve the precision and extend the 
region for form factor measurements. 



FIG. 17: The 77* — >■ |7i) transition form factor multiplied by 
SQVS in comparison with the 77* — >■ 7r° transition form fac- 
tor [l| . The dashed line indicates the asymptotic limit for the 
tt" form factor. The dotted, dash-dotted, and solid curves 
show predictions of Ref. [2y| for the asymptotic DA [30l |. 
the Chernyak-Zhitnitsky tt" DA [H, and the 7v° DA from 
Ref. [s^l, respectively. 
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FIG. 18: The 77* — > |s) transition form factor multiplied by 
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